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ABSTRACT 

Context. Direct measurement of oxygen abundance for metal-rich galaxies from electron temperature is difficult or impossible since 
temperature- sensitive auroral lines generally become too weak to be measured. 

Aims. We aim to derive the electron temperature (T e ) in the gas of metal-rich star-forming galaxies, which can be obtained 
from their ratios of auroral lines [O ii]/U7320,7330 to nebular lines [O n]/l3727, in order to establish a more robust mass- 
metallicity relationship, and compare the T e -based (O/H) abundances with those from empirical strong-line calibrations, such as 
#23 (=([0 n]i3727+[0 in]^i4959,5007)/H/3). 

Methods. We obtained 27 spectra by stacking the spectra of several hundred (even several thousand) star-forming galaxies selected 
from the SDSS-DR4 in each of the 27 stellar mass bins from log(M*)~ 8.0 to 10.6 (in units of log(M Q )). This "stack" method 
sufficiently improves the signal-to-noise ratio of the auroral lines [O ii]/U7320,7330. Using a two-zone model for the temperature 
structure, we derive the electron temperature t 2 in the low ionization region from the [O n]/Li7320,7330/[O ii]/13727 ratio, and then 
use a relation derived by fitting H n region models to estimate the electron temperature t 3 in the high ionization region from t 2 . Then, 
the direct (O/H) abundances are obtained from t 2 , t 3 and the related line-ratios. The emission lines have been carefully corrected for 
dust extinction using the B aimer line ratio after correcting for the underlying stellar absorption. 

Results. Combining our results with those from the literature with lower metallicities, we are able to provide a new relationship 
between the direct measurements of (O/H) and R23, which still shows an upper and a lower branch with the transition around 
12+log(0/H) ~8.0-8.2. It also shows that the empirical R23 method will overestimate log(0/H) by 0.2 to 0.6 dex. The new metal-mass 
relationship of the galaxies with moderate metallicities is fitted by a linear fit (12+log(0/H) = 6.223+0.23 lxlog(M*)) confirming 
that empirical methods significantly overestimate (O/H). We also derived their (N/O) abundance ratios on the basis of the T e method, 
which are consistent with the combination of the primary and secondary components of nitrogen. 

Conclusions. This study provides for the first time a method to calibrate direct O/H abundances (from T e ) for a large range of galaxies 
within a stellar mass range of ~5 10 8 M o to 4 10 10 M o . 
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1. Introduction 

The chemical properties of stars and gas within a galaxy pro- 
vide both a fossil record of its star formation history and infor- 
mation on its present evolutionary status. Accurate abundance 
measurements for the ionized gas in galaxies require the deter- 
mination of electron temperature (T e ) in the gas, which is usu- 
ally obtained from the ratio of auroral to nebular line intensi- 
ties, such as [O ra]^4959, 5007/[O m]^4363. This is generally 
known as the "direct T e -method" since T e is directly inferred 
from observed line ratios. However, it is well known that this 
procedure is difficult to carry out for metal-rich galaxies since, 
as the metallicity increases, the electron temperature decreases 
(as the cooling is via metal lines), and the auroral lines even- 
tually become too faint to measure. Instead, other strong neb- 
ular line ratios have to be used to estimate the oxygen abun- 
dances of metal-rich galaxies (12+log(0/H)>8.5), such as the 
#23 (=([0 ii]^3727+[0 iii]/U4959,5007)/H#) parameter (Pagel 
et al. 1979; Tremonti et al. 2004, hereafter T04, and the ref- 
erences therein), N2 index (=log([N n]/16583/Ha)), and some 
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other metallicity- sensitive strong-line ratios (Kewley & Dopita 
2002; Pettini & Pagel 2004; Stasinska 2006; Liang et al. 2006; 
Yin et al. 2007a). Pilyugin (2001a,b) and Pilyugin & Thuan 
(2005) also suggested their P-method to estimate the (O/H) 
abundances from the R23 and P (=[0 m]/([0 11] +[0 m])) param- 
eters. 

However, these "indirect" calibrations may produce some 
errors in metallicity estimates. For example, the double- valued 
problem of R23 for oxygen abundances, with the transition oc- 
curring near 12+log(0/H)~8.4 (log# 2 3~0.8) (McGaugh 1991). 
Most of the other strong-line ratios are sensitive to the photoion- 
ization parameter in the ionized gas (Kewley & Dopita 2002), 
and some of them suffer from the same double- valued problem 
as #23. 

Researchers continuously try to obtain the "direct" T e - 
based oxygen abundances for the metal-rich galactic gas, though 
there are some disagreements (Stasinska 2005). For example, 
Kennicutt et al. (2003; hereafter KBG03), Bresolin et al. (2004, 
2005), Bresolin (2007) and Garnett et al. (2004a,b) obtained 
the oxygen abundances from T e for a handful of H 11 regions 
in M101 and M51 with 12+log(0/H) > 8.4 by using high 
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quality spectra taken from the Very Large Telescope (VLT), 
Multiple Mirror Telescope (MMT) or Keck, 4- 8m telescopes. 
Pilyugin et al. (2006) derived a flux-flux(ff)-relation to estimate 
the [O m]/l4363/[0 m] ,15007 ratio from the obtained R 23 by 
using a sample of 48 H n regions with moderate metallicities. 
However, all these impressive progresses are only based on indi- 
vidual H ii regions. It is very difficult to measure the temperature- 
sensitive auroral lines from the integrated spectra of metal-rich 
galaxies, hence derive their "direct" metallicities. 

The wonderful Sloan Digital Sky Survey (SDSS) database 
allows this to be realized. It provides high quality spectra 
covering 3 800-9200 A with moderate resolution (~3A) for 
a large number of galaxies at different given oxygen abun- 
dances or stellar masses. It makes it possible for us to mea- 
sure the temperature- sensitive auroral lines [O n]/U7320,7330 
from the integrated light of metal-rich galaxies, which are usu- 
ally the strongest auroral lines in the measured optical spec- 
tra of the metal-rich galactic gas. To do so, we combine the 
spectra of several hundred (even several thousand) galaxies 
in each of the different mass bins (see Sect. |2] for details). 
This "stack" method significantly improve the signal-to-noise 
(S/N) ratio of the spectra, which makes it possible to use 
[O n]/U7320,7330/[O n]/l3727 ratios to derive their electron 
temperature T e in the low ionization region fe), then to obtain 
their "direct" (O/H) abundances. This is probably the first result 
concerning the derivation of "direct" (O/H) abundances from T e 
from the integrated light of metal-rich star-forming galaxies. 

This paper is organized as follows. The sample selection cri- 
teria and spectral combination are described in Sect. 2. The cor- 
rections for stellar absorption and flux measurements are de- 
scribed in Sect. 3. In Sect. 4, we re-calibrate the relations of 
the (O/H) abundances versus stellar masses and R23 parame- 
ters for star-forming galaxies on the basis of the SDSS-DR4 
database. The determinations of the oxygen abundances from 
T e are presented in Sect. 5. In Sect. 6, we compare the T e - 
based (O/H) abundances with those derived from other empir- 
ical strong-line methods, such as R23, P and N2, and derive a 
new relationship between (0/H)j; and R23. In Sect.7, we obtain 
the log(N/0) abundance ratios of these combined galaxies. The 
discussions are given in Sect. 8. We conclude the paper in Sect. 
9. Throughout this paper, the stellar masses of galaxies M* are 
in units of solar mass M Q . [O n]/L17320,7330 can also be given 
as one single line [O n] 273 25, so both of these two descriptions 
will appear in the text. 

2. Sample selection and spectral combination 

We aim to measure the electron temperatures, hence 
the 7>based (O/H) abundances, for massive galaxies. 
[O n]/Li7320,7330 lines are usually the strongest auroral 
lines in the measured optical spectra of star-forming galaxies 
(Bresolin et al. 2005), thus we focus on the T e derived from 
the line-ratio of [O n]/Li7320,7330/[O n]A3727. To increase 
the S/N ratio of the spectra, we combine the spectra of several 
hundred (even several thousand) galaxies selected following the 
criteria described below, which makes it possible to measure 
[O n]/L17320,7330 auroral lines for the massive galaxies up to 
log(Af*)=10.6. 

Firstly, we select -70,000 star-forming galaxies with 
8.0<log(M*)<10.6 from SDSS-DR4 by following their stellar 
masses and metallicities estimated by the MPA/JHU groupQ 
(T04; Kauffmann et al. 2003; Brinchmann et al. 2004; also 
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see Liang et al. (2006) for the criteria for star-forming galax- 
ies). This criterion allows us to select all the star-forming 
galaxies with available estimates of stellar masses and metal- 
licities (above the required S/N ratios) from the SDSS-DR4, 
but with stellar masses smaller than 4 10 10 M o . This mass 
limit is caused by the difficulty in obtaining measurements of 
[O n]/l/i7320,7330 auroral lines for the most massive galaxies. 

Secondly, we divide these galaxies into 27 bins with different 
stellar masses from log(M*)=8 to 10.6. The mass bins and the 
total numbers of the sample galaxies in the bins are given in 
Col. 3 and Col. 5 in Table [T] respectively. 

Thirdly, in each bin, we select about half of the galaxies 
with stronger emission lines as working samples, which allows 
us to obtain [O n]/U7320,7330 lines at a high S/N ratio for the 
galaxies with stellar masses up to log(M*)=10.6. Namely, we 
select the half of the galaxies having equivalent width values of 
[O n]/l3727, W(0 11), larger than the mean value of the whole 
sample of galaxies in each bin for the case of log(M*)<10.0, and 
larger than 30A for the case of log(M*)>10.0. The W(0 11) of the 
selected galaxies and the mean W(0 n) mean (Col. 4 in Table [T] 
given as the actual EW([0 11])) in each bin were gathered or 
calculated from the measurements published by the MPA/JHU 
group. The numbers of the selected samples in the bins are given 
in Col. 6 in Table [T] The selected sample for this study finally 
consists of 23,608 galaxies. In the following (see Sect. @]), we 
will show that the above criterion creates almost no bias in esti- 
mating the mass-metallicity relationship. 

Fourthly, the hundred (even several thousand) spectra of 
the sample galaxies in each bin are combined to be one sin- 
gle spectrum, which has high S/N ratio (-60 at continuum) 
and sufficiently strong [O ii]/LT7320,7330 lines for measuring. 
Figure 1 shows all the combined spectra at rest- frame in 27 mass 
bins, and the small figure in each panel shows the auroral lines 
[O n]/L17320,7330 region. The basic data of the combined spec- 
tra are given in Table [TJ 

These combined spectra may provide a good representation 
of the integrated light of the real galaxies in a given stellar mass 
range. Obviously, they gather more light from more individual 
H 11 regions of the galaxies with similar properties, e.g., the sim- 
ilar stellar masses. However, one may wonder whether stacked 
spectra may keep the information that underlines the real physics 
of each of the individual galaxies. Indeed the intrinsic property 
of this is similar to what happens when studying the global prop- 
erties of an individual galaxy: it is a mix, mainly of the H 11 
regions of the galaxy. Here we simply sum up a much larger 
number of H 11 regions in a large number of galaxies, and as- 
sume that the resulting properties are near the average properties 
of each individual galaxy. In Sect. 0] we will verify that this as- 
sumption indeed applies for the empirical determination of O/H 
(through the so-called R23 method), and the stacking method will 
be discussed further in Sect(Hl2. 

From our integrated-light spectra, it is still almost im- 
possible to measure other temperature- sensitive emission-lines 
[N n]/l5755, [S n]M)72. And the SDSS spectrum does not 
cover [S m]/L19069,9532 though it is possible to measure the 
[S m]/16312 emission-line for the sample with log(M*)<9.6. 
Thus it is impossible to obtain directly the r([Sm]) from the 
[S m] ratio and the T([Nn]) from the [N 11] ratio for these sample 
galaxies. In a two-zone model of temperature structure, T(N + ) is 
equal to T(0 + ) and it is acceptable to omit T(S +2 ) since it rep- 
resents an intermediate-ionization zone. 
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Table 1. The basic properties of the 27 combined spectra. Col. (1)- 
(6) refer to the counting number, the label for the mass bin (marked in 
Fig. 1), the mass bin, the average EW([0 n]/13727) of the sample galax- 
ies in each bin, the total number of the galaxies, and the number of the 
selected spectra to be combined (about half) in each bin, respectively. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


No. 


Label 


Mass bin 


EW([0 ii]) 


Total 


Selected 






(logM*) 


average (A) 


Number 


Number 



149 
158 
228 
341 
533 
593 
766 
1029 
552 
627 
705 
751 
814 
837 
906 
951 
1010 
1020 
1025 
1048 
1076 
1120 
2145 
1813 
1417 
1148 
846 



Note (a) : for the massive galaxies with log(M*)>10.10, we select those 
with stronger [O n] emission-lines with EW([0 n])< -30A for spectral 
combination. 

3. Correction for stellar absorption and extinction 

The B aimer emission lines should be corrected carefully by the 
underlying stellar absorption before measuring their fluxes prop- 
erly. This correction dominates the main uncertainties for the 
flux measurements of the B aimer lines in this study since the 
S/N ratios of the combined spectra are quite high. 

3.1. Stellar absorption correction 

The methodology used here is fitting an observed spectrum, in- 
cluding its continuum and absorption lines, with a linear com- 
bination of simple stellar populations (SSP) generated using the 
population synthesis code of Bruzual & Chariot (2003, BC03), 
GAL AXE V. The BC03 models incorporate an empirical spectral 
library with a wavelength coverage (3200-9500A) and spectral 
resolution (~3A) that is well matched to that of the SDSS data. 
We use 10 templates of SSP of BC03 with different ages (0.005, 
0.025, 0.1, 0.2, 0.6, 0.9, 1.4, 2.5, 5, 10 Gyr) at solar metallic- 
ity that are obtained by using the stellar evolution tracks of the 
"Padova 1994 library", and the Chabrier's IMF (Chabrier 2003). 
These 10 templates are a sub- sample of the 39 templates used 
by T04 to fit the spectra of the SDSS galaxies, which consist of 
three metallicities Z = 1/5, 1, 2.5Z©. The 10 templates with Z 
used in this study should be sufficient to correct the stellar ab- 
sorption for the spectra, since the degeneracy of age and metal- 
licity still cannot be perfectly broken. 

We use the program, STARLIGHT^, developed by Cid 
Fernandes and colleagues, to synthesize the stellar absorptions 
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Fig. 2. The fitting on the continuum and stellar absorption for 
one example spectrum m881. The top panel is for the spectral 
range from 3670 to 7505 A, and the bottom panel is for the de- 
tailed fitting from Ca H K lines to UJ3 (from 3910 to 4940 A). 
The solid lines indicate the observed spectrum, and the dotted 
lines give the fittings by using the SSP of BC03 and the program 
STARLIGHT (see Sect. l3.lb . The top and bottom plots are made 
by using the same data and model. 

and continua of the sample galaxies (Cid Fernandes et al. 2005, 
2007; Mateus et al. 2006, Asari et al. 2007). It is a program to fit 
an observed spectrum Ox with a model M A that adds up TV* (we 
use 10 as mentioned above) spectral components from a pre- 
defined set of base spectra. Both dust extinction on the SSP tem- 
plates and line-of-sight stellar motions modeled by a Gaussian 
distribution with dispersion cr have been considered in the pro- 
gram. The Galactic extinction law of Cardelli et al. (1989, CCM) 
with Ry=3. 1 is adopted. The fit is carried out with the Metropolis 
scheme, which searches for the minimum^ 2 . We show an exam- 
ple in Fig.[2]for the spectral fittings. 

3.2. Flux measurements and extinction correction 

We use SPLOT in the IRAlQ package to measure the fluxes of 
the emission-lines from the stellar- absorption subtracted spec- 
tra. Dust extinction (Ay, at V-band), is estimated by using the 
B aimer-line ratios Ha/H/3, Hy/H/3 and Ha/Hy, and by assuming 
case B recombination with a density of 100 cm" 3 , a temperature 
of 10 4 K, and the predicted intrinsic Ha/Uj3 ratio of 2.86 and 
Hy/H/3 ratio of 0.466 (Osterbrock 1989). Since the stellar ab- 
sorptions underlying B aimer emission-lines have been corrected 
carefully by using the synthesized SSP, the dust extinction val- 
ues derived from the three line ratios should be almost identical. 
Figure [3] shows the excellent consistency of the derived values 
for the dust extinction. 

The uncertainties of the flux measurements are dominated 
by the underlying stellar absorption correction. Another source 

3 IRAF is distributed by the National Optical Astronomical 
Observatories, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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339 
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m881 
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-71.53 
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m891 


8.9-9.0 


-65.51 


797 
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m901 


9.0-9.1 


-62.24 


1227 


6 


m911 


9.1-9.2 


-60.42 


1470 
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m921 


9.2-9.3 


-57.76 


1885 


8 


m931 


9.3-9.4 


-55.02 


2428 
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m9405 


9.40-9.45 


-52.80 


1329 


10 


m94505 


9.45-9.50 


-51.31 


1478 


11 


m9505 


9.50-9.55 


-50.41 


1670 


12 


m95505 


9.55-9.60 


-48.14 


1746 
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9.60-9.65 


-46.58 


1920 
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9.65-9.70 


-45.11 


2001 
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9.70-9.75 
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2552 


20 
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2605 


21 
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10.10-10.20 


-29.12 fl 
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24 
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10.20-10.30 


-26.28° 


6004 


25 


ml031 


10.30-10.40 


-23.79° 


6084 


26 


ml041 


10.40-10.50 


-22.18° 


5973 


27 


ml051 


10.50-10.60 


-20.72° 


5341 
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Fig. 3. Comparisons between the extinction values derived from 
Hof/H/3, Hy/HyS and Ua/Uy for the sample galaxies. A v (Ua/Uj3) 
is adopted as the actual extinction in this study. 

comes from the independent measurements by two us us, YCL 
and SYY. Figure [3] shows that the error-bars of Ay(Hy/H/3) are 
moderately larger than the other two from Ha/Hj3 and Ha/Hy, 
simply because the Hy line is weaker than the other two lines. 
Ha is the strongest of these three lines, and is affected less than 
the other two by the stellar absorption. In the following, we adopt 
A v (Ha/Hj3) as the actual extinction inside the galaxies. 

4. Consistency checks of the empirical O/H 
estimates using the R 2 3 method 

It is necessary to re-derive the relationship between (O/H) abun- 
dances and stellar masses and R23 parameters for nearby star- 
forming galaxies on the basis of the SDSS-DR4 database by fol- 
lowing the method of T04 since our sample galaxies are selected 
from DR4. This revises somewhat the calibrations obtained from 
the DR2 database by T04, but should be more reasonable for this 
study. 



Figure |4^,b shows the relationship between 12+log(0/H) 
and log(M*) and log(7?23) for the DR4 star-forming galaxies 
(-85,000, the small points). The dashed lines are the calibrations 
of T04 obtained from DR2, and the solid lines are the re-derived 
relations from DR4, given as the following formulas: 

12 + log(0/H) = -5.292 + 2.568(logM*) - 0.1 1441 (logM*) 2 , (1) 

12 + log(0/H) = 9.412 -2.251(logR23) + 4.885(logR 23 ) 2 

- 4.564(logR 23 ) 3 . (2) 

They are obtained by fitting the median values in the bins of 
0.1 dex in log(M*) and 0.05 dex in log(0/H), respectively. 

Figure |4^,b reveals the differences between the results from 
DR4 and DR2 for these two relationships which might be due 
to the expansion of the database. Eq. (Q]) will result in a lower 
log(0/H) abundance than the previous calibration at the given 
stellar mass, generally about 0.1 dex. Eq. ([2]) will result in a 
higher and lower log(0/H) abundance than the previous one at 
the given log(7?23) at the metal-rich and metal-poor ends, respec- 
tively. 

Figure[4]: shows our selected galaxies (the small blue points) 
that we extracted from the whole sample of SDSS-DR4 galaxies 
following the criteria described in Sect. [2 The large solid circles 
refer to the median-value points in the 27 mass bins given in 
Table [TJ The solid line is the same as in Fig. |4^ and Eq. (Q]). 
It shows that these median-value points of the selected samples 
are very close to the median- value relation of the whole sample 
(the solid line). Notice however that it slightly biases O/H values 
towards lower metallicities at the highest mass end. 

Figure \5\ presents the 12+log(0/H)# 23 versus log(M*) rela- 
tionship of our 27 stacked spectra. Their 12+log(0/H)# 23 abun- 
dances are estimated from the extinction-corrected R23 param- 
eters following Eq. (|2]). The corresponding stellar mass in each 
bin is the middle value there (see Table Q]). Notice that the data 
points obtained from the stacked spectra are very similar to the 
large filled circles in Fig. [4]:, i.e., the actual median- value points 
of the sample galaxies in the same mass bins. Moreover, the 
stacked spectra provide values quite close to the solid line, which 
represents the median- value relation of the whole sample. This 
implies that for the empirical O/H estimate using the R23 method, 
our stacking method provides a good representation of the me- 
dian value for the whole sample, and that the stacking method is 
able to give a reasonable estimate of the behavior of the global 
mass-metallicity relationship. 

5. The "direct" oxygen abundances derived from 
electron temperature T e 

A two-zone model for the temperature structure within the galac- 
tic gas was adopted. In this model, T e ([0 11]) (in K, as ti in 10 4 K) 
represents the temperature for low-ionization species such as + 
(and N + ), while T e ([0 m]) (in K, as t$ in 10 4 K) represents the 
temperature for high-ionization species such as ++ . 

We estimate t 2 from the ratio of [O n]^7320,7330 
to [O n]/l3727 by using the task temden. nebular in the 
IRAF/STSDAS package (de Robertis et al. 1987; Shaw & 
Dufour 1995); in the high-ionization zone is determined from 
t2 following an equation derived by fitting H 11 region models. 
Several versions of the t^-h relationship have been proposed 
(Pilyugin et al. 2006). We use the one given by Garnett (1992), 
which is valid over the range 2000K<r,([O m])<18,000K, and 
has been widely used: 

t 2 = QJt 3 + 0.3. (3) 
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i°g(v..) 

Fig. 4. Relation of the log(0/H) abundances versus the stel- 
lar masses (a) and versus the log(7?23) (b) of the nearby star- 
forming galaxies selected from SDSS-DR4 (the small points, 
-85,000). The two solid lines in (a) and (b) refer to the new 
calibrations derived from the SDSS-DR4 star-forming galax- 
ies with 8.5<log(M*)<11.3 (Eq. ©), and 8.3<12+log(0/H)<9.2 
(Eq. (|2])), respectively. The dashed lines refer to the calibrations 
obtained by T04 from DR2 database, (c) represents the sample 
of galaxies selected for this study. The large filled circles refer 
to the median-value points in the mass bins following Table [TJ 
which are very close to the solid line, the same one as in Fig.|4^. 
(Please see the online color version of the plots for the lines.) 
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8 9 10 11 

log(M.) 

Fig. 5. Relationship between the 12+log(0/H)# 23 abundances 
and the stellar masses for our 27 combined galaxies (the points). 
The solid and dashed lines are the same as in Fig.H^. 

We can then obtain the "direct" oxygen abundances from 
electron temperature for these combined galaxies by using the 
equations published by Izotov et al. (2006) for the determination 
of the oxygen abundances in H n regions for a five-level atom. 
They used the atomic data from the references listed in Stasinska 
(2005). The formulas are 

12 + \0g(0 + /H + ) = l0g(/ [0 //]i3726+i3729/%) + 

1 676 

5.961 + 0.401og^ 2 - 0.034^ + 

ti 

log(l + 1.35x 2 ), (4) 

where x^ - I0~ 4 n e t~ l/2 , and n e is the electron density in cm -3 , 
and 

12 + \0g(0 ++ /H + ) = log(/ [O ///]^4959 + i5007/%) + 

6.200 + -^-^ - 0.551og^ 3 - 0.014*3. (5) 

h 

The total oxygen abundances are derived from the following 
equation: 



The electron densities in the ionized gas of the galaxies are 
calculated from the line ratios [S n]^6717/[S n]^6731 (Cols. 2,3 
in Table 2) by using the five-level statistical equilibrium model 
in the task temden.nebular in the IRAF/STSDAS package at 
r e =10,000K. Notice that x^ in Eq. © has a very small impact 
since it is generally less than 0.1 with n e < 10 3 cm -3 . 

Table 2 lists the derived properties of the combined galax- 
ies, including the mass-bin, the [S ii]/16717/[S ii]/16731 ratios, 
the electron density, the [O n]/U7320,7330/[O n]A3727 ratios, 
the electron temperature T e ([0 n]) (in K), T e ([0 m]) (in K), and 
the 7>based 12+log(0/H) abundances, which are approximately 
8.0 to 9.0. 

Figure [6] shows the relationship between the T e -based 
log(0/H) abundances and the stellar masses log(M*) for the 
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Fig. 6. Relationship between the (O/H) abundances derived from 
T e and the stellar masses of the combined galaxies. The long- 
dashed line is the linear least- square fit for the data points. The 
solid and short-dashed lines are the same as in Fig. [5j 



combined galaxies (the filled circles). The thick long-dashed line 
is the linear least- square fit for the data point and is given by 

12 + log(0/H) = 6.223 + 0.231 x (logM*), (7) 

with an rms about 0.12dex, which shows the increasing 
log(0/H) following the stellar masses. 

This shows that the T e -based log(0/H) abundances of the 
upper branch are generally about 0.2-0.4 dex (down to 0.6 dex) 
lower than the R23 -based abundances (the solid and dashed lines, 
the same as in Fig. [5]) at a given stellar mass. The error-bars 
marked here come from the flux measurements. We did not con- 
sider the uncertainties from the conversion relationship between 
t2 and £3 , and the temperature fluctuations and gradients (see the 
discussions in Sect. [8]). 

6. Comparisons between the 7>based (O/H) and 
those from empirical strong-line methods 

In this section, we compare the T e -based O/H abundances of the 
sample galaxies with those derived from other "strong-line" ra- 
tios, such as R23, P-method and N2 (=log([N n]/Har)) index. The 
corresponding parameters and the derived (O/H) abundances are 
presented in Table 2. 

6.1. TheR 2 3 method 

The R23 parameter is widely used to estimate metallicities for 
metal-rich galaxies. The empirical relationship between O/H and 
R23 has been suggested and extensively discussed in the litera- 
ture (see Pagel et al. 1979; McGaugh 1991; Kobulnicky et al. 
1999; Tremonti et al. 2004 and references therein; Sect. [4] in 
this study). However, some recent observations show that the R23 
will overestimate the log(0/H) abundances by 0.2-0.5 dex on the 
basis of about a hundred H 11 regions in the spiral galaxies M101 
and M51 (KBG03; Bresolin et al. 2004, 2005; Bresolin 2007; 
Garnett et al. 2004a). 



Figure [7^ shows how these combined SDSS galaxies (the 
filled circles) are distributed in the relationship between log(7? 2 3) 
and 12+log(0/H)r e . The empirical R23 calibrations obtained 
by T04 (the dot-dashed line), K99 (the short-dashed lines), 
Kobulnicky & Kewley (2004, KK04; the dotted line) and that 
derived from the SDSS-DR4 database in this study (the long- 
dashed line, see Sect. |4]) are also given. The H 11 regions studied 
by KBG03 and Castellanos et al. (2002, hereafter CDT02), and 
metal-poor galaxies studied by Yin et al. (2007a) are also plot- 
ted. We will not present here other H 11 regions from Garnett et 
al. (2004a) and Bresolin (2007) since these are special H 11 re- 
gions in the center of the spiral galaxies and are not exactly the 
same case as galaxies with their integrated light of many H 11 
regions studied here. 

The combination of all these observational data with T e - 
based (O/H) estimates shows that the R23 parameters result in 
double- valued abundances with a transition region of 7.9-8.2 of 
12+log(0/H) and 0.9-1.0 of log(7? 23 ). We consistently find that 
most of our sample galaxies lie in the upper-branch with their 
8.3<12+log(O/H)<9.0and 0.4<log(# 23 )<0.75. Our samples are 
similar to those of CDT02, and more metal-rich than those of 
KBG03 and Yin et al. (2007a). Moreover, it is clear that gener- 
ally the empirical calibrations for upper-branch result in higher 
log(0/H) abundances at a given R23, up to 0.6 dex. 

We obtain a third-order polynomial fit for the total sample of 
these observational data (702, without considering the two most 
metal-rich H 11 regions from CDT02 and the three scattered ones 
with log(#2 3 )>1.0 from the SDSS) for log(7? 23 ) vs. log(0/H), 
which is given by the solid line in Fig.lTJ), and by 

log/? 23 = 73.13 - 32.78x + 4.806x 2 - 0.230x 3 , (8) 

where x=12+log(0/H). It reveals the discrepancy between this 
newly derived calibration based on the 7>(0/H) abundances and 
those derived from the previous strong-line empirical calibra- 
tions (other lines). 

To show this discrepancy more clearly and directly, we com- 
pare the (0/H)* 23 and (0/H) Te in Fig. ED In Fig.[8jt, the (0/H)* 23 
abundances of all the data points are estimated by using Eq. (|2]), 
which is the new calibration derived from the DR4 database in 
this study. In Fig. [8J3, the (0/H)# 23 abundances are estimated 
from T04's calibration derived from DR2, which could be rep- 
resentative of other empirical calibrations. They show that the 
R23 parameter will overestimate the actual O/H abundances by a 
factor up to 0.6 dex for the moderate metal-rich galaxies. For the 
objects within the transition region of abundances, the calibra- 
tions of R23 derived from DR4 will give consistent abundances 
as T e , but others (e.g. DR2) will still overestimate the O/H abun- 
dances (see the open circles in Fig.[8^,b). 

6.2. TheP-method 

It is worth while comparing the T e -based abundances with those 
derived from the P-method, which is often assumed to give "re- 
liable" metallicities for the metal-rich galaxies (Liang & Yin 
2007). Pilyugin (2000, 2001a,b; as P01) suggested his P-method 
to estimate the oxygen abundances of galaxies in which the oxy- 
gen abundance can be derived from two parameters, R23 and P 
(=[0 m]4959,5007/([O n]3727+[0 m]4959,5007)). He derived 
the 0/H=/(7?23, P) formulas from a sample of metal-poor H 11 
regions with 7.1<12+log(0/H)r e <7.95 and a sample of moder- 
ately metal-rich H 11 regions with 8.2<12+log(0/H)r e <8.7. The 
fitting relation for moderately metal-rich H 11 regions (~40) is 
given as Eq. (8) in Pilyugin (2001a). Pilyugin & Thuran (2005; 
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Fig. 7. (a) The relation of the T e -based (O/H) abundances with 
R23 parameters for our sample galaxies (the filled circles). 
Correlation with those from the literature: the large open cir- 
cles refer to the H 11 regions in M101 taken from KBG03, the 
squares represent the H 11 regions given by CDT02, and the 
small open circles refer the low-metallicity H 11 regions and 
galaxies taken from Yin et al. (2007a). Some empirical cal- 
ibrations are also given: those of T04 (the dot-dashed line), 
Kobulnicky et al. (1999; K99, the two dashed-lines, y=log0 32 = 
log([0 m]4959,5007/[O n]3727)), Kobulnicky & Kewley (2004; 
KK04, the dotted line), and that derived from SDSS-DR4 in this 
study (DR4, the long-dashed line), (b) The third-order polyno- 
mial fit (the solid line) for all the observational data points given 
in (a), also given as Eq. (8]). Other calibrations from the literature 
are the same as in (a). 



as PT05) have renewed the P calibrations by including sev- 
eral improvements, such as enlarging the sample (-104 metal- 
rich H 11 regions). However, on the basis of -20,000 metal-rich 
star-forming galaxies from the SDSS-DR2, we found that the 
oxygen abundances derived from this revised P-method (PT05) 
are ~0.19dex lower than those derived from the previous one 



Fig. 8. Comparison between the T e -based and the R23 -based 
log(0/H) abundances of the sample galaxies: (a) the R23- 
calibration was derived from DR4 with 12+log(0/H)>8.3 in this 
study (Eq. (|2])); and (b) the R23 -calibration was derived from 
DR2 by T04 with 12+log(0/H)>8.5. The symbols as the same 
as in Fig. [7] The solid line in each panel represents the equal re- 
lation, and the two dashed lines are the relations by increasing 
12+log(0/H)# 23 abundances by 0.2 and 0.6 dex, respectively. 



(P01), and in particular are -0.60 dex lower than the Bayesian 
abundances obtained by the MPA/JHU group (Yin et al. 2007b). 
These abundance discrepancies strongly correlate with the P pa- 
rameter and weakly depend on the log(7?23) parameter, which is 
especially serious for the case of P < 0.55 (Yin et al. 2007b; 
PT05, P01). 

In Fig. [9j we compare the 7>based log(0/H) abundances of 
the 27 combined galaxies with those derived by using the P- 
method of P01 for the metal-rich branch. The H 11 regions stud- 
ied by KBG03 and CDT02 are also plotted. It seems that the 
P-method will almost result in a constant (O/H) abundance of 
about 12+log(0/H)~8.4 for both of our samples and most of the 
data from the literature, although the log(0/H)/> and log(0/H)j> 



8 



Liang et al.: Oxygen abundances of metal-rich galaxies from T e 




8 8.5 

12+log(0/H)(Te) 



Fig. 9. Comparison between the T e -based (O/H) abundances and 
those estimated from P-method for our sample galaxies. The 
symbols as the same as in Fig. [71 
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Fig. 10. Comparisons between the T e -based (O/H) abundances 
of the sample galaxies with those estimated from N2 index cal- 
ibrations: (a) from D02; (b) from PP04; (c) from Y07; (d) the 
(N2,P) calibration from Y07. 



are consistent within 0.3 dex. The reason for our sample galax- 
ies may be that the POl's P-method is not quite reliable to de- 
rive their oxygen abundances. This is because all our moder- 
ate metal-rich sample galaxies selected from the SDSS have 
P <0.55 (see Table 2), due to the [O m] line becoming weaker 
in the metal-rich environments, which is the case for most of the 
metal-rich SDSS galaxies (see Yin et al. 2007b). But the calibra- 
tion formula of P01 (as well PT05) was derived from fitting the 
sample H n regions that mostly have P >0.55, and only a few of 
their samples have P <0.55. The samples of CDT02 have simi- 
lar P- values to ours, and show similar trends to those described 
here. The samples of KBG03 generally have higher P values 
than 0.55. The reason for their discrepancy from the equal-line 
may be that they are just in the turn-over region of abundances 
when strong-line calibration cannot give accurate abundances. 

6.3. The N2 method 

The [N n]/Ha emission-line ratio (as N2=log([N n]/l6583/Haf)) 
is useful to estimate the metallicities of galaxies (Liang et al. 
2006; Yin et al. 2007a, as Y07; Pettini & Pagel 2004, as PP04; 
Denicolo et al. 2002, as D02; Kewley & Dopita 2002), though it 
may depend on the N-enrichment history of the galaxies (Liang 
& Yin 2007), and ionization parameters (Kewley & Dopita 
2002). N2 is not greatly affected by dust extinction due to the 
close wavelength positions of [N n] and Ha, and the devel- 
oped near infrared spectroscopic instruments can gather these 
two lines for the galaxies with intermediate and high redshifts, 
which stand in the early epoch of the universe. 

Here we compare these T e -based (O/H) abundances with 
those derived from the N2 calibrations of D02, PP04, Y07, also 
the (N2,P) calibration from Y07 (P =[0 m]/([0 n]+[0 in])), 
given in Fig.[T0h-d. They show that the ^-calibration given by 
PP04 could give (O/H) abundances that are more consistent with 
the T e -based ones, while other ^-calibrations often overesti- 
mate the (O/H) abundances for these moderate metal-rich galax- 
ies. 



7. The log(N/0) abundance ratios 

The "direct" electron temperature in the low-ionization region 
of the galactic gas and the high quality optical spectra make it 
possible to estimate the log(N/0) abundances of these combined 
galaxies properly, which is useful to understand the "primary" 
and/or "secondary" origin of the nitrogen element. 

If the "seed" oxygen and carbon for the production of nitro- 
gen via the CNO processing are those incorporated into a star at 
its formation and a constant mass fraction is processed, then the 
amount of nitrogen produced is proportional to the initial heavy- 
element abundance, and the nitrogen synthesis is said to be "sec- 
ondary". If the oxygen and carbon are produced in the star prior 
to the CNO cycling (e.g. by helium burning in a core, followed 
by CNO cycling of this material mixed into a hydrogen-burning 
shell), then the amount of nitrogen produced may be fairly inde- 
pendent of the initial heavy-element abundance of the star, and 
the synthesis is said to be "primary" (Vila-Costas & Edmunds 
1993). In general, primary nitrogen production is independent 
of metallicity, while secondary production is a linear function of 
it. 

We adopt the low-ionization region temperature t<i for 
[N n] emission regions fe([N n])= ^([O n])) and the 
[N n]/U6548,6583/[0 n] ,13727 ratio to derive the log(N/0) 
abundances of these combined galaxies by using the formulas of 
Izotov et al. (2006) for the log(^) and log(^r) abundances, with 

§ = §r- Figure [Upresents the results. These 7>based (O/H) and 
(N/O) abundances show that the sources of nitrogen for these 
moderate metal-rich galaxies are consistent with the combina- 
tion of the primary and secondary components. This is a robust 
investigation on the production of nitrogen for these moderate 
metal-rich galaxies since both the N and O abundances are ob- 
tained from electron temperature on the basis of their integrated 
light. The H n regions of KBG03 and CDT02 also show a similar 
nitrogen source to our samples. Most of the dwarf galaxies taken 
from van Zee & Haynes (2006) show the primary source of their 
nitrogen. 
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Fig. 11. Log(N/0) abundance ratios of the sample galaxies as 
functions of their 12+log(0/H) abundances derived from T e (the 
large filled circles). The dot-dashed line refers to the "primary" 
component of nitrogen, the long-dashed line refers to the "sec- 
ondary" component of nitrogen, and the solid line refers to the 
combination of these two components, which are taken from 
Vila-Costas & Edmunds (1993). Other recent observations are 
also given: the triangles represent the dwarf galaxies studied by 
van Zee & Haynes (2006), and the open circles and squares 
represent the extragalactic H n regions studied by KBG03 and 
CDT02, respectively. 



8. Discussion 

Our procedure to derive T e -based oxygen abundances for the 
metal-rich galaxies may be affected by some biases that need 
to be discussed further. The major contributors for such biases 
are the dust estimates, the stacking method, and the uncertain- 
ties related to the temperature determination. 

8.1. Dust extinction 

Our method is strongly dependent on the dust extinction, since 
the determination of the temperature depends on the ratio of red 
over blue emission lines. However, the S/N of the stacked spec- 
tra is so large that our estimate of Ay is very accurate. On the 
other hand, one may wonder whether by stacking spectra, we 
still present the representative of the mean properties of individ- 
ual galaxies. 

We confirm the strong correlation between the dust extinc- 
tion (Ay) and the stellar masses of star-forming galaxies, which, 
as expected, means the more massive galaxies generally have 
more heavy dust attenuation (see also Salim et al. 2007). In 
Fig. [I2l the large filled circles represent our stacked galaxies, 
the small points represent the whole sample of the star-forming 
galaxies from SDSS DR4 (same as in Fig.|4^), and the large open 
circles represent the median values of this whole sample in the 
corresponding mass bins. In most of the considered stellar mass 
range, the stacked spectra provide a Ay-M* relationship almost 
identical to that of the whole sample, meaning that the stacking 
method and selection criteria do not cause significant biases in 
the estimation of the dust extinction. Notice however that at the 



Fig. 12. Relationship between dust extinction and stellar masses 
of the sample galaxies. The large filled circles refer to the results 
from the stacked spectra of the galaxies in the 27 mass bins. The 
small points refer to the star-forming galaxies from SDSS DR4 
(the same as in Fig.|4^), and the open circles refer to the median 
values of them. (Please see the online color version of the plot 
for the points, especially for the open circles.) 



highest mass range, stacking spectra provide a slightly higher ex- 
tinction, which may well be due to the selection method. Indeed, 
by selecting stronger star forming galaxies than the average, it is 
not surprising to find slightly higher dust extinction. 

8.2. Are stacked spectra providing a good representation of 
rea! gaiaxies ? 

The "stacking" method has been used widely in astronomical 
study, e.g., Zheng et al. (2006) and Chang et al. (2006), who 
stacked the images or property parameters of galaxies to study 
their properties. There are also examples of uses of stacked spec- 
tra: Mathis et al. (2006) stacked the spectra in several mass bins 
to study the cosmic star formation history; Baldry et al. (2002) 
analyzed the stacked low-resolution spectra of 166,000 galax- 
ies at redshift 0.03<z<0.25 extracted from the 2dFGS; Schiavon 
et al. (2006) used the stacked Keck DEIMOS spectra to study 
the ages and metallicities of some red galaxies; and Wild et al. 
(2007) used the stacked SDSS spectra in redshift bins to study 
the [O ii]/13727 nebular emission from the galaxies hosting Ca n 
absorption line systems and galaxies hosting Mg n-selected DLA 
absorbers. 

We are however conscious that the stacking method may 
produce undesirable biases since it is unclear how some phys- 
ical properties may be affected by this method. In Sect. |4] and 
Sect. [SI 1 , we have convincingly shown that the stacking method 
does not cause significant biases in representing the median 
value for the O/H estimated values from the R23 method, and for 
the extinction estimates, respectively. Both estimates are based 
on ratios of lines at different wavelengths. It is thus tempting 
to believe that temperatures, estimated in a similar way, would 
not be affected by the stacking method. However, we agree that 
this needs a full demonstration, which is beyond the scope of 
this paper. We also notice that the integrated spectrum of a sin- 
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gle galaxy comes from all the contributed star-forming regions. 
By using stacked spectra from different galaxies, we are study- 
ing the resulted light coming from a larger set of star-forming 
regions. In principle our method is not less accurate than the de- 
termination of the O/H through temperature for a single galaxy, 
and thus may not alter the calibrations proposed in this paper. 

The major interest of the adopted stacked method presented 
here, is that it allows us for the first time to sample the stel- 
lar mass versus log(0/H) relation using a temperature dependent 
method, up to relatively high stellar mass galaxies. 

8.3. The uncertainties of the T([0 ll]A7325)-based 
abundances 

In the optical band, electron temperatures could usually 
be derived using five sets of auroral/nebular line inten- 
sity ratios: [O m]^4363/[0 m]^4959, 5007, [S m]^6312/ 
[S m]^9069, 9532, [S u]U63 17,673 1/[S n]^4072, 
[N n]^5755/[N n]^6548,6583, and [O n]^7320,7330/ 
[O n]/l3727. However, for the integrated spectra of the SDSS 
sample galaxies with moderate stellar masses, only the in- 
dicator [O n],U7320,7330/[O n]^3727 is available to derive 
their electron temperatures in the galactic gas since only 
[O ii]/U7320,7330 lines (marked as A1325 sometimes) are 
strong enough to be measured from these spectra. Therefore, 
the line ratio [O n]/U7320,7330/[O n]A3727 provides the only 
"direct" measurements for the electron temperature and the 
"direct" oxygen abundances for these moderate metal-rich 
galaxies, though some arguments exist to doubt the accuracy 
discussed as following. 

One of the problems of Til 325) is the contribution of re- 
combination to the [O ii]/17325 line. However, this has been es- 
timated by KBG03 using the formula given by Liu et al. (2000), 
who found that recombination typically contributes less than 5% 
to the [O n] line flux, which corresponds to a temperature er- 
ror of only ~2%-3% or less than 400K for their worst case. A 
correct treatment of recombination should take into account the 
effect of temperature gradients within ionized nebulae. 

The second problem is the inconsistency between the 
7X7325) and the 7(5755) and 7(4072), presented by Bresolin 
et al. (2005). By using the VLT/FORS spectra, Bresolin et al. 
(2005) obtained the 7(7325) for 32 H n regions located in some 
spiral galaxies, and 7(5755) and 7(6312) for only half of their 
sample since these two lines are much weaker. In addition, 
[S n]4072, from which 7(4072) was derived, has been seldom 
detected, since its measurement is made difficult by low S/N 
ratio in the spectra. They compared these temperature values, 
and found that 7(7325) seems to overestimate the temperature if 
compared to 7(5755), while the opposite happens for 7(4072). 
However, the error-bars of the latter two temperatures are much 
larger than 7(7325), and these two lines are much weaker than 
A1325. We cannot measure the other temperature- sensitive auro- 
ral lines from most of our sample spectra. 

Another problem is that the 7 [On] vs. 7[0m] relation 
does not match the model fitting relation very well (see Fig. 1 of 
KBG03; Garnett 1992). They mentioned that the sources of this 
disagreement remain unresolved, but one of the sources may be 
the recombination contribution (Liu et al. 2000), others could be 
the radiative transfer effects and the observational uncertainties. 
However, their objects with large scatter mostly have 7[0 m] be- 
tween 9000K and 12000K, and there may be a correlation be- 
tween 7 [On] and 7[0m] in the extremes of the temperature 
range, e.g. T e <9000K and T e >13000K. Most of our sample 
galaxies have just 7[Om]<8000K, thus no direct observations 



are against that our galaxies follow that correlation of 7[0 n] vs. 
7[0m]. 

Therefore, the temperature derived from [O \\]A1325 is ac- 
ceptable, and is the only way to get the "direct" oxygen abun- 
dances for metal-rich galaxies. 

8.4. Temperature fiuctuations and temperature gradients 

It has been argued for many years that directly measured elec- 
tron temperatures and the corresponding abundances from col- 
lisionally excited lines may have systematic errors due to tem- 
perature fluctuations and temperature gradients, especially for 
the metal-rich H n regions (Peimbert 1967; Stasinska 1980; 
Stasinska 2005). However, the auroral lines have been used to 
derive the "direct" electron temperature and the oxygen abun- 
dances for the metal-rich H n regions in recent years (KBG03; 
Garnett et al. 2004a,b; Bresolin et al. 2004,2005; Bresolin 2007). 
KBG03 has discussed very clearly the "reliability" of the auro- 
ral line method for abundance estimates. We simply summarize 
some arguments and discussions presented by KBG03, Garnett 
(1992) and Stasinska (2005) regarding this topic. 

In the presence of fluctuations in T e , forbidden line strengths 
can overestimate the average temperature, and therefore under- 
estimate the true nebular abundances. However, the question of 
the existence of temperature fluctuations is still controversial 
(Peimbert 1967; Dinerstein et al. 1985; Shaver et al. 1983). The 
large-scale temperature gradient mimics the effects of temper- 
ature fluctuations, and can over-estimates the real temperature, 
and consequently the oxygen abundance is systematically under- 
estimated (Garnett 1992; Stasinska 2005). 

Despite these concerns, as KBG03 discussed, there are good 
reasons to believe that the forbidden lines give abundances that 
are close to correct, at least for the range of O/H of interest 
here. Recent measurements of infrared forbidden lines (that have 
weak temperature dependences) in the spectra of planetary neb- 
ulae yield abundances in good agreement with those from op- 
tical forbidden lines, even where there is a large discrepancy 
with recombination line abundances (Liu et al. 2000, 2001). 
Likewise, measurements of radio recombination line tempera- 
tures of Galactic H n regions give results that are consistent with 
forbidden line measurements of the same objects (Shaver et al. 
1983; Deharveng et al. 2000). In both cases, the line emissivity 
is only moderately dependent on T e . It would be worthwhile to 
explore the metal recombination lines and infrared fine- structure 
lines for the determination of electron temperature in the extra- 
galactic H ii regions of high metal content. However, there are 
obvious difficulties in exploiting these techniques: recombina- 
tion lines are very weak and hard to detect in extragalactic H n 
regions (Esteban 2002), and infrared observations require orbit- 
ing telescopes (Garnett et al. 2004b). 



9. Conclusion 

We have derived direct oxygen abundances based on electron 
temperature for metal-rich galaxies, using an original method 
of stacking. To do so, we select a large sample of star-forming 
galaxies (-23,608) from the SDSS-DR4 and obtain 27 spec- 
tra by combining the spectra of the several hundred (even sev- 
eral thousand) galaxies within 27 various stellar mass bins from 
log(M*)=8.0 to 10.6. This stacking method improves the S/N ra- 
tios in a great deal of the spectra, without providing major biases. 
These high quality data with moderate resolution (3A) make it 
possible for us to measure the [O n]/Lr7320,7330 emission-lines 
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for the galaxies with stellar masses up to log(M*)~10.6, then de- 
rive their electron temperature in the low ionization region by 
the ratio of [O n]/U7320,7330/[O n]A3727, hence the T e -based 
O/H abundances. 

It provides a new calibration of the 7?23 method, which 
may well apply to galaxies with metallicities ranging from 
12+log(O/H)~7.0 to 9.0. It also shows that the empirical R 2 3 
calibrations overestimate the log(0/H) abundances by about 
0.2 to 0.6 dex. We also derive a new relationship between T e - 
derived (O/H) abundances and stellar masses from these moder- 
ate massive galaxies, which can be fitted by a linear fit, namely 
12+log(0/H) =6.223+0.231xlog(M,). The 7>based log(N/0) 
abundance ratios show that their nitrogen sources are consistent 
with the combination of the primary and secondary components 
for these moderate metal-rich galaxies. 
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Fig. 1. The 27 combined spectra of the sample galaxies at rest-frame in the corresponding mass bins. The small figure inside each 
panel shows the auroral lines [O n]/U7320,7330. The label for the mass bin of each spectrum is given in the bottom right corner. 
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Table 2. The derived properties of the combined galaxies. 
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Notes: Col. (1)-(17) consequently refer to: (1) the mass of the center value of the bin; (2) the ratio of [S n]/l6717/[S ii]/16731; (3) the electron density; (4) the ratio of 
[O n]/U7320, 7330/[O n]/l3727; (5) the election temperature in low ionization region; (6) the election temperature in high ionization region; (7) the oxygen abundances 12+log(0/H) 
derived from election temprature; (8) the dust extinction A v estimated from Ha/Hfi; (9) logarithm of the ratio of ([O ii]/13727+[0 iii]/U4959, 5007)/H/3; (10) logarithm of the ratio of 
[O m]AA4959, 5007/[O n]i3727; (11) the 12+log(0/H) abundances obtained by using the R 2 3 calibrations derived from DR4 (Eq.[2j; (12) the 12+log(0/H) abundances obtained by using 
the R 2 3 calibrations derived from DR2 (T04); (13) the P parameter [O m]^4959, 5007)/ ([O n]^3727+[0 m]^4959, 5007); (14) the 12+log(0/H) abundances obtained by using the 
P-method calibration of P01; (15) the N2 index (=log([N n]/l6583/Hor)); (16) the 12+log(0/H) abundances obtained by using the N2 calibration of PP04; (17) the log(N/0) abundance 
ratios derived from election temperature method. 



